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Eifective medium theory of transport in disordered systems, whose basis is the replacement of spatial disorder 
by temporal memory, is extended in several practical directions. Restricting attention to a 1-dimensional system 
with bond disorder for specificity, a transformation procedure is developed to deduce, from given distribution 
functions characterizing the system disorder, explicit expressions for the memory functions. It is shown how 
to use the memory functions in the Lapace domain forms in which they first appear, and in the time domain 
forms which are obtained via numerical inversion algorithms, to address time evolution of the system beyond 
the asymptotic domain of large times normally treated. An analytic but approximate procedure is provided 
to obtain the memories, in addition to the inversion algorithm. Good agreement of effective medium theory 
predictions with numerically computed exact results is found for all time ranges for the distributions used except 
near the percolation limit as expected. The use of ensemble averages is studied for normal as well as correlation 
observables. The effect of size on effective medium theory is explored and it is shown that, even in the asymptotic 
limit, finite size corrections develop to the well known harmonic mean prescription for finding the effective rate. 
A percolation threshold is shown to arise even in 1-d for finite (but not infinite) systems at a concentration 
of broken bonds related to the system size. Spatially long range transfer rates are shown to emerge naturally 
as a consequence of the replacement of spatial disorder by temporal memories, in spite of the fact that the 
original rates possess nearest neighbor character. Pausing time distributions in continuous time random walks 
corresponding to the effective medium memories are calculated. 

PACS numbers: 5.60Cd, 61.43.-j,72.80.Ng 



I. MOTIVATION FOR THE STUDY 

Description of the movement of excitations and quasiparti- 
cles is crucial to the study of a variety of disciplines in physics 
and allied sciences Conductivity in metals and semi- 
conductors, energy transport in molecular aggregates, atomic 
motion in ceramic materials, molecular hopping in cell mem- 
branes, all present a diverse variety of contexts in which such 
description is indispensable. At a sufficiently macroscopic 
level, the description is often provided by a Master equation 
of the type 

dPJt) 



dt 



(1) 



which governs the evolution of the probabilities P,„(f) of oc- 
cupation of site m by the moving entity at time t via the tran- 
sition rates F„,„. Here m is typically a vector index in the 
appropriately dimensioned space. The method of solution 
of such an equation relies on the diagonalization of the so- 
called A-matrix defined through A,„„ - -F„i„ for m + n and 
Amm - En Fnm- One Can always formally write the solution of 
the probability vector P{t) from its initial value P{0) as 

P{t) = e-^'PiO). (2) 

If the system is ordered (quasiparticle moving on a crystal lat- 
tice), this solution becomes practical because the diagonal- 
ization can be performed via discrete Fourier transforms and 
the A;th mode of the probabilities f *(f) = 2™ PmiOe''"" can be 
written as 



Inversion into explicit P,„(f)'s is straightforward, and the 
specifics of the dynamics of the A-matrix are evident through 
the eigenvalues A'' of the latter. An alternative way of an- 
alyzing the transport is via random walks. The relationship 
between a random walk description and a Master equation de- 
scription was given long ago by Bedeaux et al. |l3|]. 

Physical systems in which disorder cannot be neglected are 
rampant in nature. Whether the lack of order arises because 
some transfer bonds are weaker or stronger than others, or 
whether the quasiparticles must surmount barriers at some lo- 
cations but not at others, disorder must be seriously taken into 
account in the description of these systems. A natural way is 
to replace the given system by a corresponding ordered prob- 
lem with temporal memory. What this means is that the orig- 
inal time-local disordered problem, given by Eq. ([T), is re- 
placed by the so-called generalized master equation (GME), 



dPJt) 
dt 



^,„„(t-s)P„(s)-'Wn,n{t-s)P,n{s), (4) 



where the memory functions 'W,„„ are of the form ^,„_„, 
i.e., translationally invariant, the replacement equation being, 
therefore, soluble via discrete Fourier transforms. Then, in 
the Laplace domain (e is the Laplace variable and tildes de- 
note Laplace transforms), the counterpart of Eq. (3) is 



P'(e) 



(5) 



p*(f) = e-^''P''(0). 



(3) 



where y[„,„ - -1Vm„ for ni + n, JKmm - li„ 'W,™,, and J^k is 
the Fourier transform of Jim-n- Generalized master equations 
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were introduced in the sixties to understand the derivation of 
the irreversible Master equation from reversible quantum me- 
chanics and a comparative review of methods has been given 
by Zwanzig [4]. Needless to say, depending on one's choice, 
one could employ, for the description of transport in disor- 
dered systems, continuous time random walk (CTRW) for- 
malisms (?] instead of GME formalisms in equations such 
as (|4]i. The two ways of description, CTRW's and GME's have 
long been shown 17[] to be entirely equivalent to each other in 
fully arbitrary (rather than space-time decoupled) forms since 
1974. The space-time decoupled case appeared in Kenkre et 
al ||6D and the demonstration of the equivalence for the general 
(arbitrarily coupled) form was given in Kenkre and Knox |7] 
(see, for instance, their Eqs. (43)). It appears that the gen- 
eralization in ref. |7] was missed by several different authors 
who reported it independently but subsequently [8„ even 
six years later ifiotl . Some of this was commented on in brief 
in ref flT] (see Eqs. (27-31) of the latter reference). 

Two questions are important to answer in the context of this 
program of the description of a disordered system. To what 
extent is the replacement of the spatial disorder by temporal 
memories possible and meaningful even in principle? And 
what is the prescription to calculate the memories and effec- 
tive transfer rates given appropriate information about the dis- 
order in the particular system? Without the first, it is senseless 
to begin. Without the second, the study is useless. 

The first question can be answered quite trivially on a little 
reflection. Consider Eq. (1) solved. By assigning the so- 
lutions for the probabilities, P,„(r), to an appropriate ordered 
lattice, carry out the direct Fourier transform to obtain P''{t). 
Put the Laplace transform of the latter into 



(6) 



and Fourier and Laplace invert to get the (translationally in- 
variant) memories J?l„,„, equivalently "Wmn appearing in the 
GME dU. The presence of initial conditions in the above 
prescription means that each possible set of initial conditions 
would have a corresponding set of memory functions, a situ- 
ation which is obviously unacceptable for practical purposes. 
However, in order to turn Eq. (|6|l into a practical prescrip- 
tion for computing memories which is independent of initial 
conditions, all that is necessary is to carry out an ensemble 
average over the possible realizations of disorder, compati- 
ble with what is known (for instance a distribution function) 
about the disorder. Such an average makes the system transla- 
tionally invariant after the average. Then the first term in the 
right hand side of Eq. (|6]l which is the reciprocal of the Fourier 
and Laplace transform of the (ensemble averaged) propagator, 
is independent of initial conditions. The propagators directly 
lead to the memories. 

This is precisely the method devised long ago by Kenkre 
|fl2ifl3|] to obtain exact expressions for memory functions an- 
alytically for a quantum mechanical (not disordered) system, 
although no ensemble average was involved in that context. 
Because equations ([T) and (|4]i as well as the operation of av- 
eraging over configurations are linear, it quite unnecessary to 
make any assumptions or offer demonstrations to be able to 



state with certainty that the replacement program is possible. 
Analyzing the problem from the viewpoint of the application 
of projection techniques ifTill to the problem, it also becomes 
clear from Zwanzig's formal theory that a memory will auto- 
matically appear in a closed description of any quantity that is 
formally projected from another whose evolution equation is 
time-local. Here the projection is represented by an ensemble- 
average over disordered realizations. This too requires no cal- 
culation, only a moment's reflection. The initial condition 
problem, rarely discussed in the disorder context, also makes 
its appearance in the projection formalism lfl4l 1511 . It appears 
as a separate term. In the original context l\M it is removed 
through the initial random phase or diagonality assumption. 
In our present disorder context it disappears on carrying out 
the ensemble average we have mentioned above. 

What is really necessary in the sense of calculations comes 
to the second question we have posed above, i.e., the find- 
ing of an explicit and practical prescription that would allow 
one to go from information about the disorder in the real sys- 
tem to the memories (or pausing time distribution functions) 
in the replacement problem. Very few instances of such a 
prescription exist in the literature, a noteworthy attempt be- 
ing in the early work of Scher and Lax [16] that gave sup- 
port to the well known theory of Scher and Montroll [17]: 
on a phenomenological basis, the latter addressed with great 
success unexplained puzzles of transport in xerographic ma- 
terials. Known information about transfer rates between ran- 
domly located sites is converted via an approximation scheme 
in the appendices of ref. (l^, into the continuous time random 
walk pausing time description. That is the kind of prescription 
that one needs in developing a usable theory. 

The present paper focuses on a different manner of con- 
verting disorder information into temporal memories that has 
to do with the venerable subject of effective medium theories 
(EMT) iillillll^lIIlllllzillTlllillilll. We provide 
the essential background on EMT in section 2, along with a 
prescription we provide in a particularly convenient form that 
transforms the disorder into explicit memory functions via a 
double Laplace transform procedure. Our prescription facil- 
itates the extraction of the new results we present in subse- 
quent sections. The spirit of the investigations we present is 
most akin to, among early attempts that have discussed mem- 
ory functions in the EMT context lfloll23ll . the work of Haus 
andKehr[^25J. 



II. EXPLICIT DISORDER-TO-MEMORY TRANSFORM 

Effective medium theories are unabashedly approximate, 
i.e., do not claim to provide an exact solution of the problem. 
They sacrifice exactness for practicability, i.e., prefer useful- 
ness to avoidance of approximations. One of the first instances 
of their application is by Bruggeman [19]j3ut many later and 
independent reports may also be found I2CII|. The basic idea, 
explained in many textbooks and reviews fW\, is to assume 
that the memory represents an effective ordered medium in a 
mean field sense with a magnitude (of the memory) which is 
such that any departures, introduced in keeping with whatever 
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information is known about the disorder, average out to zero, 
thus estabhshing the ordered system as representing a varia- 
tionally optimum Hmit. 

To understand this concretely, we consider from now on in 
this paper a 1 -dimensional case of Eq. ^ with bond disorder. 



dPrM 

dt 



f„,+ l[P,„+l(f) - P,„(0] + F„[P„,-i(t) - P,n{t)-\, (7) 



the disorder being expressed via a distribution p(/). What 
this means is that any transition rate F„, can have any pos- 
itive value / with probability density p(f) normalized such 
that p{f)df - 1 . No correlations exist in the actualization 
of rates at different locations. The replacement of the disor- 
dered time-local system by an ordered system with memory 
then proceeds by writing in place of Eq. Q, a GME 



dt 



f 

Jo 



dsT(t - s)[P,n^i(s) + P„,_i(5) - 2Pm(s)], (8) 



which is translationally invariant and describes elemental 
transfer interactions that are nearest neighbor as in the orig- 
inal (disordered) problem. 

When applied to the present system, the general effective 
medium concept requires the following procedure. One eval- 
uates the probability propagators for two different systems: 
the ordered system obeying Eq. ([8]l, and a system obeying 
Eq. dHJ augmented by terms that represent a single disor- 
dered bond. Transport across that bond occurs not through the 
memory Tit) but through a rate / drawn from the distribution 
pif). One averages the latter propagators over the distribu- 
tion, i.e., carries out an integration of the result with weight 
pif), and demands that the average equal the corresponding 
ordered propagators. Details may be found in the references 
given and lead straightforwardly to 



df 

Jo 1 + 2[/ ■ 



pif) 



r(e)][^o(e)-^i(e)] 



= 1. 



(9) 



The above equation is an implicit equation for the memory 
T that can be obtained in principle from the given probability 
distribution function pif), the quantities if/o and being the 
propagators of the ordered system: the probability of remain- 
ing on the site initially occupied is if/Q whereas the probability 
of occupation of the adjacent site is i/'i. Equation® is the 
same as Eq. (22) of ref. iH, or Eq. (5^ of ref. or Eq. 
(7) ofref. [23], orEq. (3.17) of ref. lUi or Eq. (38) of ref. 
rt22ll . Through a simple manipulation we rewrite it first as 



Jo 



df 



Pif) 



^ [l - e«Ao(e)] (10) 



and then, by introducing a quantity f , in the remarkably sim- 
ple and convenient form 



Jo f + € f + { 



(11) 



The quantity ^(e, T) is a function of both e and of 'Fie) 
since the selfpropagator (^o depends explicitly on 'F(e) as well 



as on e. Generally, 



^ = Tie) 



eipoie) 



1 - e-Ao(e) 



(12) 




For the infinite 1-d chain with nearest neighbor rates, given 
that, in this case, etfio equals [1 + 4;F(e)/e]"'^^, one has the 
specific expression 



(13) 



As we will see below, this restatement ( fTTT i of the basic 
EMT equation (|9]) allows us to obtain a number of our re- 
sults in a straightforward fashion. With very few exceptions 
in the literature, the result (|9|l is used in the long-time limit 
and therefore involves the Markoffian replacement of Fit) by 
^(0[Jq Fis)ds]. This is equivalent to the e — > limit. By an 
Abelian theorem eijfo becomes identical to the f ^ oo limit of 
the selfpropagator which is zero if the system considered is an 
infinite chain. Equation ( fTTT i then reduces to the well known 
effective medium theory result Ii29il that the effective transfer 
rate Fe/f - Fis)ds = Tie 0) equals the harmonic 
mean of the disordered /'s: 



1 _ 1 _ r 
''~f " ¥iO) " Jo 



df 



pif) 
f ■ 



(14) 



By contrast, our interest in the present paper is to extract 
new information from the memory equation without taking 
the Markoffian limit, and to go beyond common uses of effec- 
tive medium theory. We will derive some general features of 
the EMT memory in section 3, describe our extensions of the 
theory for times that are not asymptotic in section 4, analyze 
the emergence of spatially long range memories in section 5, 
study finite size effects in section 6, and present conclusions 
in section 7. 

The first of the results of our present investigation is the re- 
formulation implicit in Eq. (fTTT i interpreted as a transform of 
the distribution function pif) (disorder information) into the 
effective medium quantity Fit) (temporal memory). Specifi- 
cally, we can regard Eq. ( fTTT i as related to a double Laplace 
transform. One applies the direct Laplace transform twice: 
first to pif), with a dummy variable y as the Laplace variable, 
to obtain giy), and then to giy) with ^ as the Laplace variable 
to obtain hi^): 



giy) 



Jo 



pif)e-yfdf; hiO 



Jo 



giy)e ^^dy. 



The prescription for extracting the memory Fit) in the 
EMT equation (|8]l from the disorder distribution pif) consists, 
thus, of computing the double transform hi^) of the disorder 
distribution, equivalently performing the integral on the left 
side of (fTTT i, and inverting into the time domain the memory 
transform Fie) after solving for it from the implicit equation 



hi^) = 



1 



Fie)+^ie,F) 



(15) 



One has, thus, a practical prescription to obtain the f depen- 
dence of the memory Tit) from the disorder p(/). 

The usefulness of the form of the basic equation we have 
presented, Eq. ( fTTT i. should be already clear by comparison 
to the well-known asymptotic result for the effective rate Eq. 
(fl4l l and noticing from Eq. (fT2l l or Eq. (fTsT l that in the asymp- 
totic limit ^ vanishes. Further uses are reported in the rest of 
the paper. 



m. NATURE OF THE EFFECTIVE MEDIUM MEMORY 
FUNCTIONS 

We have applied the prescription of Eq. ( fTTT i to various 
distribution functions p{ f) to obtain Tie) and discovered that 
the results share a number of common features. These features 
become apparent on inverting the transform to obtain Tit), the 
memory function in the time domain, and can be understood 
as we show below. The numerical inversion scheme we use is 



n(f,M) 



In 2 

t 



^ -/ ln2\ 



minftM) :M+\lM\ll;\l 1 \ 



where f^git, M) is the approximation to the Laplace inverse of 
Tie). The precision required to sum the series, i.e. the num- 
ber of significant digits, is 2.2M while the precision of the 
resulting expression is 0.90M. So, if one uses double preci- 
sion in the calculations, the value of M should not be larger 
than 7. For a detailed discussion on the numerical inversion 
of Laplace transforms, see ref. |34]. Note that the function is 
only evaluated at the real and positive values of the Laplace 
variable e. 

The result of the Laplace inversion is always that the EMT 
memory Tit) consists of two pieces, a 5-function at the ori- 
gin of time (f = 0) and a part that is negative but finite. A 
schematic depiction is given in Fig. 1. In order to understand 
this and other qualitative, and some quantitative, aspects of 
Tit) from general arguments, consider the actual system evo- 
lution equation ^ on the one hand, and the representative 
EMT equation ([8]) on the other, both for an initial occupation 
of only the site m. Let us first evaluate the first time derivative 
of Pm{t) at the initial time. The respective results are 



dPJt) 



dt 



- -(Fm + Fm+l) 



(17) 



(=0 



for the actual Master equation, and 



dPM 



dt 



^0+ 

= I d 

/=() -'o- 



sTit - s)[P,n+l(s) + Pm-lis) - 2P,„(s)] 

(18) 

for the representative EMT equation. A configuration average 
over the distribution p(/) converts the right hand side of Eq. 
dnli into -2</) = -2 / dfp(f)f. It is impossible for Eq. (US 




time 

FIG. 1 : Shape of the effective medium memory function 'F(t) show- 
ing the delta function of strength {/> at the origin and the negative 
piece Q(t). The time integral from to oo of the memory TU) is 
l/<l//>. 



to yield a non-zero result (because of the limits of integration) 
unless f{t) contains a ^-function at the origin. As Eqs. ( fTTl ) 
and ( fTSl l must yield results that equal each other, we deduce 
that the form of the EMT memory function is 



nt) = imt) - Q(t). 



(19) 



The origin of the 5-function at f = is clear from the above 
analysis. That the additional part must have a time integral for 
all time which is negative follows from the general result ( fT4l i 
that the integral over all time of fit) is the harmonic mean 
1 1 {I If) - J dfp{f)/f which is always smaller lUsfl than the 
arithmetic mean (/). 

We also note that the integral of Q(t) over all time is now 
determined: 



f 

Jo 



Q(t)dt = </) 



1 



(20) 



JO <1//) 

Additional information can be obtained in this exact manner 
about the memory function, for instance, the initial value of 
Q(t). Differentiation of Eq. (|7]i with respect to time yields the 
initial second time derivative 

]d^P,M] ^ 2 (^-2^^ + Fl + F,n^iFj) . (21) 



dt^ 



1=0 



Similarly, differentiation of the EMT generalized master equa- 
tion (ISj yields, after a configuration average. 



dt^ 



= 6{fy + 22(0). 



(22) 



(=0 



Carrying out the configuration average of the former result, 
which gives 4(/^) + 2(/)^, and equating the two values of the 
second time derivatives at the initial time, we can evaluate 
2(0) exactly for any distribution function as 



2(0) = 2[{f}-{ff] 



= 2 



(/.«/)/=)-(/./.(/./) 



(23) 
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It is also straightforward to continue in this manner with fur- 
ther differentiations to obtain exact initial values of higher 
derivatives of Q{t). For instance, in terms of the A-matrix 
appearing in Eq. (|2]i, we can evaluate the initial value of the 
rth derivative of P,„ via 



d'Pn 



df 



and proceed as shown above with configuration averages. 
For our nearest neighbor rate system we have Amn - 

We do not carry out this program here but use the limited 
information gathered above to develop a simple analytical ap- 
proximation to the memory cast in the form of a difference 
of a term proportional to a delta function and another to an 
exponential. This 'exponential' approximation to the EMT 
memory for any given distribution of the rates is 



Tait) ^ {f)6{t) - 2[{f) - {ffyy 



(24) 



The subscript a clarifies that the memory is approximate. 
While the precise shape of the actual memory function will 
not be captured by our approximation (l24b . examples to be 
given in the next section will make clear that the approximate 
memory can be remarkably good. 

The general behavior of the time dependence of the mem- 
ory function consisting of a decay (infinitely fast for our sys- 
tem) to negative values and then a rise which is often rela- 
tively slower is typical in many systems. It is usually encoun- 
tered in studies of the velocity autocorrelation (v(f)v) which is, 
needless to say, very closely related (in our case simply pro- 
portional) to the memory function. The small time behavior 
represents initial transfer at a higher rate; the subsequent be- 
havior is affected by disorder or imperfections in the system 
as they are encountered in the motion. Indeed, the velocity 
autocorrelation for a random walker completely confined to a 
finite space exhibits this very behavior, the overall integral of 
{v{t)v) for all time being precisely zero because of the con- 
finement: the mean square displacement saturates in this case 
(see, e.g., ref. ll36ll for a nuclear magnetic resonance context.) 



IV. BEYOND ASYMPTOTIC DESCRIPTION 

Research reported on the basis of effective medium the- 
ory is almost entirely focused on an asymptotic description 
of quantities such as the diffusion constant. There have been 
a few exceptions as in the work with focus on the ac con- 
ductivity, e.g., by Odagaki and Lax [26.1 , and others [18 1, the 
publications of Haus and Kehr ll23ll24l l25ll . the anisotropic 
disordered systems studied by Parris ||27'1, and the granular 
material stress work of Kenkre [30]. One of the questions we 
address in the present paper is how well effective medium the- 
ory works for times outside of the long and short time asymp- 
totic domains, i.e., for all times. To address this problem we 
now compute the memories explicitly for three different dis- 
order distribution functions p(/), use those memories in the 



GME to calculate observables and compare the results to nu- 
merically obtained exact solutions of the disordered Master 
equation. 



A. Some Specific Distributions p(/) 

A natural distribution to consider is the multi-delta distri- 
bution 

M 



p(f) = J]ai6{f-fd 



;=1 



wherein the nearest-neighbor transition rates may take one of 
M values f, each with a weight a,, with Yjfii ffi = 1. We will 
focus particularly on the case M = 2, so that 



pif) = a6(f - fi) + (I - a)6if - f2). 



(25) 



The arithmetic and harmonic means of the rates, additionally 
the mean of the square of the rates, are given by 

</> - a/i +(l-a)/2, 
1 /l/2 



<1//) a/2 + (l-«)/r 
(f) = aff + (l-a)fl 



(26) 



The distribution p{f) for this case is shown as the two arrows 
in Fig. 2. 

The second particular distribution we consider is the 
gamma distribution (related closely to the Poisson distribu- 
tion): 



P(f) 



— f"e-^^. 



r(n + ly 



(27) 



The arithmetic and harmonic means of the rates, and the mean 
of the square are 



<1//) 
</') 



n + 1 



r 

(«+ l)(« + 2) 

? ■ 



(28) 
(29) 
(30) 



A plot of p(/) itself is displayed for the particular case of n = 
1 and 7 = 4 in Fig. 2. 

The third case we consider is the triangular distribution 
given by 



p(/) = 



if-fo + ft)/fi 



fo-fh<f<fo- 



i-f + h + fb)lfl fo<f<h + fb, (31) 

elsewhere. 



The minimum possible rate is /o - fb and the maximum pos- 
sible rate is /o -i- fh. The distribution rises linearly from the 
minimum value with slope 1 /f^ until it attains the value I /ft, 



6 



gamma 




FIG. 2: Examples of probability distributions p(/). The two arrows 
represent tlie double-delta distribution with equal weight a = l-a = 
0.5. A gamma distribution with n = 1 and 7 = 4 and a triangular 
distribution with /o = 0.3 and b = 0.2 are depicted by the solid 
and dashed lines respectively. Units of / in the plot are arbitrary 
and the same as those of /o and b for the triangular distribution, and 
reciprocal to those of y for the gamma distribution. 



at / = /() and then descends with the same magnitude of the 
slope down to the maximum value. The meaning of /o is that 
it is the value of / at the apex (and hence the mean of the dis- 
tribution), and fb is half the length of the base of the triangle. 
The distribution is shown in Fig. 2 for /o = 0.3 and fi, = 0.2. 
It leads to 



if) = /o, 



if) 



<1//) 



(32) 
(33) 
(34) 



B. Evaluation of Memory Functions for Specific Cases of p(/) 

The approximation to the memory given by the formula 
(l24l i is easily evaluated for the three distributions by substi- 
tuting in the formula the respective values of (/), (/^) and 
1/(1//). As one example, note that for the gamma distribu- 
tion it is given by 



Tait) = 



(35) 



The memory function 'f'(t), whether derived from ( fTTT i 
or the simpler (l24l i. can be used immediately to calculate 
other, more directly observable, quantities. A useful quan- 
tity is the (dimensionless) mean square displacement (m^) - 
^mtrP'Pmit) for initial localization at the origin. It is simply 
twice the double time integral of the memory: 

{m^}^yyPM = 2 f ds f ny)dy. 



The time-dependent diffusion coefficient D(t), a quantity often 
used in transport theory to describe the instantaneous state of 
motion, may be defined as one half the product of the square 
of the intersite distance a and the time derivative of the mean 
square displacement. It is proportional to a single time integral 
of the memory function: 



2\dt 



hi: 



T(s)ds. 



These have exact expressions in terms of (/) and Q{t) appear- 
ing in Eq. ( fT9] l. If we use our simple exponential approxima- 
tion for Q(f), they become 



(</>-(i//r'f 



2 
D{f) 



[l-e-'!'), (36) 
(37) 



{II, f) 2«/2)-(/)2) 

= </)-(</)-(i//r')(i-^"'^1' 

where the time constant t is given by 

^ (/)-«i//)r' 

2«/2) - </)2) 

It is straightforward to get expressions particular to the distri- 
bution functions chosen. As expected, the mean square dis- 
placement starts out linearly with slope twice the arithmetic 
mean of the rates and ends up also linearly with slope twice 
the harmonic mean of the rates. Correspondingly, the time 
dependent diffusion constant decays from a higher to a lower 
value. 

There are a number of ways the above simple analysis can 
be put to use to extract physical information. For instance, 
the mean square displacement of a walker initially localized 
at a single site will first grow linearly but then saturate to a 
finite value at long times if there are broken bonds in the 1-d 
infinite system. Broken bonds correspond to a p(f) that has a 
non-zero value at / = which means that there are bonds at 
which the transition rate is zero. In such a case, 1/(1//), the 
harmonic mean of the rates, and consequently the long time 
D{t), vanish. Equation ( |36] | can then be used to extract the 
value at which the mean square displacement saturates at long 
times: 



Iim(OT ) 



{fr 



if) - {fv- 



(38) 



This consequence of the exponential approximation ( l24b to 
the memory is simply a case of the general EMT result 



hm(m ) 



-2 lim 



dQ{e) 



e^o de 



(39) 



This may be proved from the Laplace transform of the non- 
delta part of the memory in Eq. ( fT9l ) via a Taylor expansion: 



(l/f) 



■ -e 
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de 
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In the presence of broken bonds in 1-d, the harmonic mean of 
/'s vanishes. Since the mean square displacement is twice the 
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FIG. 3: Comparison of the exact (numerical) memory function with 
the approximate EMT given by our formula i24\ for a double-delta 
distribution in which the concentration of the smaller of the rates is 
0.9 and the ratio of the transition rates is 10. Plotted are the ex- 
act memory function (solid line) and the exponential approximation 
(dotted line). It is surprising how close the agreement is, given the 
coarse nature of the approximation. 

double time integral of 'T'(t), the limit e ^ and the use of an 
Abelian theorem establish Eq. ( [39] l quite generally. If 'F(t) is 
expressed via the exponential approximation ( l24l i. the general 
result reduces to Eq. ( l38b . 

Despite what appears as an impressive agreement of the ex- 
ponential approximation that we see displayed in Fig. [3] for 
a double delta distribution with a - 0.9, and /2//1 = 10, the 
approximation generally will not capture the actual decay in 
time for all distribution functions and may be regarded only 
as a highly simplified manner of description. For greater ac- 
curacy than can be provided by the relatively coarse approxi- 
mation of Eq. ( l24b , it is necessary to return to the prescription 
of Eq. ( fTTT i, calculate !F(e) through the solution of the implicit 
equation, and then invert the transform to obtain the memory. 
When the EMT memory is calculated in the Laplace domain 
via our prescription based on Eq. ( fTTT ), the derived quantities 
D{t) and (m^) can be obtained very simply in the Laplace do- 
main by dividing !F(e) by e and (except for proportionality 
constants) respectively. 

The calculation of T' from Eq. ( fTTI ) is easy and analytically 
doable for the double-delta distribution. Defining the quantity 

'? = (!- a)fi + a/2, 
we get the soluble cubic 

P-2r^{2T^^/e + fi+f2)) 

+r{Wif2/£ + (2/1/2 + (/ + fif - r]^)) 

-4/V|/e - (2/i/2(/ + /2) - 277/1/2) = 0. (40) 

Standard analytic formulae lead to the appropriate solution 
which can then be numerically Laplace-inverted. We have 
carried out such a procedure in the next section. 



FIG. 4: Comparison of EMT predictions with exact results for three 
disorder distributions p(/). Plotted is the time-dependent diffusion 
coefficient (normalized to its initial value) as a function of the di- 
mensionless time t = (/) t for three cases of p{ f) and a localized 
initial condition. Solid lines are the effective medium results and 
open circles correspond to the numerically exact results obtained by 
averaging over 20000 different realizations of the disordered chain. 
Dashed lines on the right show the asymptotic values of the diffusion 
coefficient. From top to bottom, p(f ) is the double-delta distribution 
with fi 7/2 = 0.5 and a = 0.5, a triangular distribtution with /o = 0.3, 
//, = 0.2, and a gamma distribution with n = l.The agreement is 
striking for all three cases, not only for extreme limits but for inter- 
mediate times as well. 



Similar procedures can be used for the gamma distribution 
and the triangular distribution. Explicit polynomials do not 
result for !F in those cases but the equations can be solved nu- 
merically and inverted. We have carried out these procedures 
for these two distributions as well and report the results after 
inversions into the time domain using Eq. ( fT6l ). 



C. Comparison of EMT and Exact Solutions for All Times 

We now display the results of the predictions of effective 
medium theory and the numerically obtained exact evolution 
not only for long times as is usually done, but for short and 
intermediate times as well. For each distribution we calcu- 
late the exact and the full EMT results. The exact results are 
obtained via numerical matrix operations as explained else- 
where in this paper With the exception of single-run studies 
to be reported further below, we repeat the operations tens of 
thousands of times, each time using a different realization of 
the chain. Then we average over all runs to produce the quan- 
tity we desire. The effective medium theory prediction for that 
quantity is also determined via the effective medium memory 
function both in its full form as given from our Eq. ( fTTI ). and 
from our analytic approximation, Eq. ( l24b . 

We first treat the case when a single site is initially fully 
occupied. 
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1. Localized Initial Condition 

Figs. |4]and|5]display the comparison graphically, the quan- 
tity selected being the time-dependent diffusion coefficient 
normalized to its initial value: D{t)/D{0). All three distribu- 
tions are considered in Fig. |4] The agreement of the effective 
medium theory (solid lines) with the exact evolution (open 
circles) is remarkably good for all cases considered and for all 
intermediate times as well. The description appears thus ex- 
cellent for the parameters considered for times that need not 
be asymptotic. 

In order to explore parameter values for which the agree- 
ment may not be as good, we restrict ourselves to the double- 
delta distribution in (a) of Fig. |5] take the two possible rates 
fi and /2 to occur with equal weight (a = 0.5), but vary the 
ratio: /1//2 = 0.5,0.1,0.01 as we go down the graph. EMT 
is still found to provide a fine description for all times but de- 
viates more as the rates become more disparate. To drive this 
situation to an extreme where the EMT would serve worst, we 
consider a broken bond system in (b) of Fig. |5] This means we 
take fi - and fi+Q for different values of the concentra- 
tion a. The large time value of D{t) is now zero and the mean 
square displacement (m^) (proportional to the integral of D(t)) 
saturates. Physically, the saturation value measures the size of 
clusters (separated by broken bonds from other clusters) on 
which the walker is localized at long times. We have already 
obtained analytic expressions for the saturation value from the 
full EMT (in Eq. ( [39] l) and from the exponential approxima- 
tion (in Eq. (l38Tl). Fig. 5b is an attempt at looking at EMT 
in the worst possible light by comparing the time evolution of 
{irp-) predicted by it to that given by exact calculations. We do 
this for the broken bond case { f\ - 0) for two concentrations a 
of broken bonds: 0.01 and 0.1 as shown. The main display in 
Fig. 5b shows the two {n-p-) curves. To make the discrepancy 
of the saturation value particularly clear, we show the inset in 
which the one case a - 0.01 is displayed on a semilogarith- 
mic scale. The abscissa is the dimensionless time {f)t as in 
the main figure. The ordinate is {m^} on a linear scale, the 
values displayed as 0.5 and 1 being 5000 and 10000 (i.e., in 
units of lO'*). The accumulated values of the mean square dis- 
placement, the localization cluster sizes, are S.VlxlO-' from 
the exact calculations but only 5.00x10^ from the EMT: both 
are denoted by solid lines in the inset. The corresponding val- 
ues for the a - 0.1 case are 88.6 and 49.5 respectively. The 
exponential approximation to the EMT is way off as it pre- 
dicts 99 for the a = 0.01 case and 9 for a - 0.1. (The latter 
is denoted by a dotted line in the inset.) This is to be expected 
from the crudeness of that approximation. 



2. Extended Initial Conditions for Single Runs 

An actual experiment in a real physical situation is per- 
formed not on an ensemble but on an individual system. How 
can EMT, which has at its root an ensemble average, pro- 
vide a valid description for the experiment? Standard Gibbs- 
Boltzmann arguments do not help as an answer here because 
our interest in using EMT is not only for asymptotic times 




FIG. 5: Worst-case scenario comparison of EMT and exact results. 
In (a) we plot the time-dependent diffusion coefficient (normalized 
to its initial value) as a function of the dimensionless time t = (/> / 
for the double-delta distribution function for a = 0.5 and, from top 
to bottom, /1//2 = 0.5, /1//2 = 0.1, /1//2 = 0.01. While good, 
the agreement gets worse for disparate /'s. To explore a regime in 
which the agreement is bad, in (b) we consider two broken-bond sys- 
tems (the ratio of the f's being zero and therefore extreme) with two 
different concentrations a = 0.1,0.01 as shown. Plotted is the mean 
square displacement showing saturation at long times. Here /i = 
and /2 = 0.2. Open circles correspond to the exact (numerical) so- 
lution obtained by averaging over 20000 different realizations of the 
disordered chain which consists of 801 sites. Solid lines are theoret- 
ical results from the EMT. The inset shows the a = 0.01 case, the 
ordinate being on a linear scale in units of 10*. See text for discus- 
sion. 



when the system might have completed the mixing process 
but for all times. One possible answer to this question might 
lie in the nature of the initial condition. If it is extended in 
space, various configurations of transition rates in a random 
system may be realized even at short times. With this idea 
in mind we now describe our investigation of extended initial 
conditions for single runs. In particular, we study the agree- 
ment of EMT and single-run evolution of the actual system as 
we vary the spatial extent of the initial condition. 

We carry out calculations from exact numerical considera- 
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FIG. 6: Relative difference between exact results and EMT, or error 
of the EMT, for a single run. Plotted as large open circles is the 
error (see text for definition) as a function of the number of sites 
initially occupied, i.e., the value 2/< + 1. No averages are performed. 
The error is seen to decrease as the initial width increases, allowing 
the walker to sample different configurations. Inset: Comparison of 
D(t)/D(0) curves for two different values of the initial width, ju = 5 
(crosses) and /i = 50 (open circles), with the EMT (solid line). It 
is clearly seen that EMT agrees with the simulations for spatially 
extended initial conditions without averaging. 



tions for systems of 801 sites without changing the configu- 
rations of the transition rates once set in accordance with the 
double-delta distribution, and take initial conditions that are 
not of the form PmiO) = 6^^, but of the extended form 

1 

PmiO) = y 6n,r 

2u + l ^ 

r=-il 

which represents a patch initial condition of spatial extent of 
2fi + 1 sites. We call this the initial width. The limit jj - 
gives us back the initial condition we have used in the studies 
above. We find that larger patches result in smaller deviations 
of the EMT predictions from the exact results. 

The inset of Fig. 6 shows values of D{t)/D{0) for a single 
configuration for two different values of the width (open cir- 
cles represent fi = 50, crosses represent ju = 5) along side the 
corresponding prediction of effective medium theory (sohd 
line). The integrated difference between the EMT result and 
the exact results depicted in that inset (as plotted on a loga- 
rithmic time scale) provides a convenient measure of the error. 
We thus define, for each value of jU, a measure of the relative 
error, through the expression 



U — oo 



(.v)-D^^(0 



L>^^ (s) 



ds 



where s = In ({f)t) . A plot of the (numerically evaluated) 
relative error as a function of the initial width /x is presented 
in the main graph in Fig. 6, and clearly shows that the relative 
error decreases monotonically as the patch width increases. 



3. Correlation Type Observables 

There are, in general, different kinds of observables that can 
be computed from the solution to the Master equation. Simple 
observables O are those which associate with each site (state) 
m, a value Om that the observable takes when the particle is in 
that state. The mean value associated with such an observable 
at any time t can then be written as 

{O it)) = Yj 0„P„ it) = Yj OA,n W Pn (0) , (41) 

m m,n 

where in the second form we have expressed the result in 
terms of the propagators *P, and the initial probability distribu- 
tion governing the particle's occupation of the possible states 
of the system. This can be put in the form 



where 



(Oit)} 



(42) 



(43) 



is the mean value of the observable given that the particle 
started in state n at time t - 0. Simple observables can thus be 
calculated by incorporating into the averaging process an av- 
erage over the different possible starting locations of the par- 
ticle. 

Correlation type observables, also of great interest in sta- 
tistical physics, do not correspond to (simple) observables of 
this type. Indeed, they span two or more different states (or the 
same state at two different times). An example is (A (f) B (0)) 
given by 

2 ^A.n it) BnPn (0) = (0>„5„/'„ (0) . (44) 

m,n n 

Here and are, respectively, the values of A and B when 
the particle is in the state m, and 



(A ii))„ — ^ Am*Pm,« it) 



(45) 



is the mean value of A at time t if the particle started in state n 
aXt -0. Consider, for instance A^ with components 



A,„ — 6n 



(46) 



It is an indicator observable taking the value 1 if the particle 
is at site { and the value otherwise. Then the correlation 
function 



{A' it) A' (0)) = (0/'r(0) 



(47) 



is just the propagator i/'^/' weighted by the relative initial prob- 
ability of finding the particle in the state {' . This shows that 
the propagators il/t^c themselves can also be considered as ob- 
servables of the system. Of course, in a specific disordered 
system, the self-propagatori^f f (f) , e.g., will depend on the lo- 
cation of site { in the disordered chain. The effective medium 
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FIG. 7: Effect of varying patch width of the initial condition. Plotted 
is the time dependence of the self propagators in a disordered chain 
(the same disordered chain is used in all of the calculations) whose 
transfer rates are drawn from a double-delta distribution with f\ jf^ = 
0.1 and a - 0.5. Self propagators at 401 sites, 200 to the left and 200 
to the right of origin are calculated. The dashed line is the average 
of all of the 401 self propagators obtained this way. The solid lines 
correspond to averaging over 3, 9, 15 and 25 (from top to bottom) of 
the self propagators around the origin. 

propagator *Po (0 may not, therefore give a good approxima- 
tion to any given self-propagator i/^f / (f) in any single realiza- 
tion of the disordered system. We intuitively expect, however, 
that self-propagators, averaged over an initial distribution of 
starting positions on the same chain, will approach that of the 
effective medium, as the width of the initial distribution of 
starting sites is increased, i.e., that 




To verify this intuition, we present calculations in Fig. 7 of 
the self-propagator averaged over such an initial distribution 
of starting sites, for different values of fi = 3,9, 15, and 25. 
A comparison of the limiting curve with the predictions of 
effective medium theory is given in Fig. 8. 

V. SPATIALLY LONG RANGE MEMORIES 

The original evolution equation that describes the disor- 
dered system, Eq. Q, is local in time and nearest-neighbor 
in the character of its transition rates. Given that effective 
medium theory provides an approximate rather than exact de- 
scription of the actual dynamics described by Eq. O, one 
may ask whether the introduction of non-locality in time in 
the EMT should be accompanied by non-locality in space as 
well. This is a natural question to pose because of the emer- 
gence of spatially long range memories that were found long 
ago lfT2l [13I1 when GME's were calculated for quantum me- 
chanical systems by a method of Eq. similar in spirit to 



FIG. 8: Self propagators obtained from the disordered chain (nu- 
merically, represented by the open circles), by the effective medium 
theory (solid line) and from an ordered chain with transfer rates F^ff 
throughout (dashed line). 

the one followed in the EMT. Stated differently, the question 
we ask is whether the replacement of Eq. ^ by Eq. ([8]) with 
nearest-neighbor transition memories is sufficient or whether 
the latter should span longer distances. The answer is pro- 
vided in this section. 

Consider Eq. (|7]i solved for ?„,(€), the Laplace transform 
of the probability of occupation of the mth site in terms of the 
matrix A'' corresponding to the configuration /u (a particular 
realization of the transition rates / throughout the system). 
Carrying out the average over the configurations ju one gets a 
translationally invariant situation: 

P"'^'^-I.{7h^) (48) 

Performing a discrete Fourier transform, we get /'*(e)/f *(0) 
which we substitute in Eq. (|6]l to get the exact memory func- 
tion: 

Because the configuration average has been carried out al- 
ready at this point, we do not display the superscript /i on the 
A. 

There is no guarantee whatsoever that the A;-dependence of 
is of the form (1 - cos k). The exact memories need not, 
therefore, have nearest neighbor character. The nature of the 
disorder will influence the A;-dependence. It is therefore clear 
that spatially long range memories will naturally develop, in 
general, their precise form being determined by the particular 
distribution p{f). 

The exact procedure is to be contrasted with the EMT pro- 
cedure, which, as explained in Section 2, necessarily results 
in the absence of spatially long range memories. This is so 
because one assumes the memories to be nearest-neighbor in 
character, and obtains them variationally. 
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FIG. 9: Spatially long range memories obtained from exact numeri- 
cal considerations plotted as a function of the Laplace variable. Units 
on both axes are of {/>. Plotted in (a) is the exact Tiie) (solid line), 
calculated for rings of 100 sites, and the almost identical EMT mem- 
ory (open circles) along with the asymptotic rate 1/(1//) (dashed 
line). Plotted in (b) on a scale blown up by almost 3 orders of mag- 
nitude are the much smaller long range memories '7^,(f) in dotted, 
dashed and solid lines, for n = 2, 3, 4, respectively. The distribution 
is double-delta with /1//2 = 0.1 and a = 0.5. 



in magnitude relative to the nearest neighbor 'Fi(e). This is 
shown clearly in (b) where the longer range memory trans- 
forms, Tnie), are shown. The scales in the plots in (a) and (b) 
differ by a little less than 3 orders of magnitude so it is indeed 
clear that the long-range memories are small. It is thus that 
the EMT can successfully describe the evolution even though 
it possesses only nearest-neighbor memories. Note that, while 
J^xie) is sigmoidal in shape, the long-range memories seem to 
peak for intermediate e and to be negligible for both large and 
small e. 



VI. FINITE SIZE EFFECTS 

To the best of our knowledge, effective medium consider- 
ations have been used only on infinitely large systems in the 
past. We present below useful EMT results for finite rings 
of sites, i.e., chains obeying periodic boundary conditions. 
The self-propagator for such a system is given in the Laplace 
domain by 



1 



2r(e)(l -cosA;) 



(50) 



where k takes on the values (27r/A^)[0, 1,2, ...A^ - 1]. In the 
long time limit, the self propagator iAo(0 tends to 1/A^ as one 
knows both from the explicit limit of Eq. dSOl l or from the 
physical statement that the probability equalizes over the ring 
sites. This means via an Abelian theorem that etfro(e) l/N 
as e — > 0. The use of this limit in Eq. (fT2b leads to an impor- 
tant long-time consequence of our general equation ( fTTT l. 



1 



N-l 



P(f) 

f + F,ff(i±jy 



(51) 



which is an extension to finite systems of the well-known 
harmonic mean result of Eq. ( fT4b . Here we have used 
Fgff = ;F(e — > 0) as earlier Equation ( BTl i must be solved for 
Feff implicitly and becomes explicit only as — » 00 when 
the Feff term within the integral disappears. 

The implicit equation for the case of the double-delta dis- 
tribution function of Eq. ( |25] ), 



In Fig. (|9]l we display the result of the full numerical exact 
procedure outlined above caiTied out on a chain of 801 sites, 
making sure during each run that the value of P,„(r) is negli- 
gible (comparable to the precision of the machine used) at the 
boundaries of the chain. The distribution used is double-delta, 
the two rates are in the ratio /i //2 = 0.1 and the concentration 
of each is equal to the other We plot in (a) the Laplace trans- 
form of the nearest-neighbor memory obtained from the ex- 
act procedure (solid line), fiie), as a function of the Laplace 
variable e, both the abscissa and the ordinate being expressed 
in units of the average rate (/). Also plotted is the result of 
the EMT procedure (dots) and the dashed line that represents 
the asymptotic rate 1/(1//). There is hardly any difference 
in the exact and the EMT result. What this must mean is 
that the non-nearest neighbor memories must be much smaller 



Feff - 



N-l 



N 



1 -a 



/l + — /2 + — 



can be converted into a quadratic equation and solved explic- 
itly. With 



J = /i(l-A? + A?a) + /2(l-A?a), 



one has 



eff 



(52) 



Normally, i.e., when both /i and /2 are non-zero, there is a 
unique solution as we discard the negative root because F^jf 
must be real. 
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FIG. 10: Bifurcation of tlie effective long time transfer rate for a 
double delta distribution in a finite system of N sites. Plotted is P^jj 
as a function of the concentration of broken bonds (i.e., bonds with 
the rate f\ = 0), the rate associated with the remaining fraction 1 - 
a of unbroken bonds being equal to /2. A transcritical bifurcation 
occurs when a equals 1 IN . For concentrations higher than this value, 
the effective rate vanishes but changes linearly with the concentration 
for lower a. Solid (dotted) lines denote the stable (unstable) solution. 



If one of the two possible rates, e.g. /i, is zero, i.e., if bro- 
ken bonds exist in the finite system, an interesting situation 
arises, both roots being of physical interest. The lower root 
is zero, not negative, in this case. If one varies the concentra- 
tion a of the broken bonds, a transcritical bifurcation occurs 
as displayed in Fig. [TO] at the point at which a equals the 
reciprocal of the number of sites in the ring. As this num- 
ber increases, the bifurcation point moves towards vanishing 
concentration. We recover the known result that, for an infi- 
nite system, the effective rate is zero for any concentration of 
broken bonds. Additionally, we get a percolation threshold for 
finite systems. The two solutions exchange stability at the crit- 
ical concentration {a - \ /N), there being transport throughout 
the ensemble-averaged system for broken bond concentrations 
below the critical value. 



'-^ 0.125 



CD 
> 




FIG. 1 1 : The behavior of the relative difference between the nearest 
neighbor effective rates calculated from the EMT and an exact nu- 
merical procedure (see text) as a function of fi //a for a sum of two 
delta functions with a = 0.5 



sites we have the exact values, lim^^o T^f^i^) 
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16/2 
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124r(l + r2) + 23Qr^ + 17(1 + r^) ' 
r(3 + 2r)(2-l-3r)(l +4r)(4H-r) 



(53) 



■ (7 + 3r)(3 + lr)(r + 1)(7 + 36r + Ir^) ' 

The effective medium memory function is given by Eq. (|52| | 
with a - 1 /2. In order to quantitatively examine how different 
the exact and effective medium values are, we define a relative 
difference as 



1 

Ji 



and plot it as a function of /i //2 in Fig. 1 1 for = 3, 4, and 
5. We find that the values predicted by the effective medium 
theory are slightly different from the exact values. The relative 
difference between the two decreases as the number of sites in 
the ring becomes larger and larger Therefore the effective 
medium theory predicts the correct values in the limit e — > oo 
when N —> CO, but finite size effects exist otherwise. Note that 
for finite A^, effective medium theory always predicts larger 
values than those that are calculated exactly. 



It is interesting to see how the effective medium nearest 
neighbor memory function compares with the exact nearest 
neighbor memory function as e ^ oo in finite rings. One 
can obtain the exact memory functions for a ring of A^ sites 
through Eq. (49) by averaging over all possible configurations 
of the ring. For simplicity, we will consider the double-delta 
distribution with a - 1/2. For rings with N - 2, 3,4, and 5 



VII. CONCLUDING REMARKS 

The purpose of this paper is to make a contribution to the 
description of transport of quasiparticles such as electrons, 
excitons, atoms, interstitials, vacancies, or other more for- 
mally regarded random walkers in a disordered system such 
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as a solid, a photosynthetic system, or a molecular aggregate. 
We have focussed our attention on incoherent motion as de- 
scribed by a Master equation and have restricted our analysis 
to nearest-neighbor transfer on a 1 -dimensional chain, infi- 
nite or finite obeying periodic boundary conditions (ring). Our 
general goal is to translate information about static spatial dis- 
order of the given system into dynamic temporal features of a 
representative ordered system, to do it explicitly by convert- 
ing distribution functions into memory functions, to study the 
validity of predictions of the effective medium theory, and to 
report several extensions we have made of the theory. 

Star ting wi th the spirit of effective medium theory used by 
many IMll IIS llltll 111126] as expressed through Eq. (Hi, 
we arrive at a transparent prescription Eq. ( fTTT i. The prescrip- 
tion is in the form of a transform, is a natural generalization 
to all times of well known results such as the harmonic mean 
recipe of Eq. (fT4l) for effective long time rates, and involves 
the solution of an implicit equation for the Laplace transform 
of the EMT memory. 

We show how to obtain the memory in the time domain by 
numerical inversions of the Laplace transform produced by 
the solution of the implicit equation. Additionally we derive, 
exactly, partial information about the memories. We provide 
an understanding of the special feature of EMT memories that 
it consists of two pieces of which one is a delta function. We 
derive a simple approximate formula in the time domain, Eq. 
( l24l i. for the memory. It can provide a rough and sometimes 
adequate representation of the exact evolution as Fig. 3 shows 
for the particular distribution and parameter set that we have 
used in that case. Lest one develop a false confidence in this 
coarse approximation, we have shown Fig. 5b in which its 
predictions for a broken bond system are quantitatively quite 
different from the exact answers. 

We use effective medium theory to go beyond an asymp- 
totic description and compare the EMT description with exact 
(numerically obtained) predictions. The quantity we choose 
for comparison is the time-dependent diffusion coefficient 
D{t) which is proportional to the time derivative of the mean 
square displacement or equivalently to the time integral of the 
memory function. Not only do we find excellent agreement at 
long and short times as expected from previous work, but sur- 
prisingly good agreement at intermediate results also. Figs. 4 
and 5 show this clearly. We carry out this comparison in two 
parts: by doing an ensemble average over initial conditions 
for localized initial placement of the walker; and by doing a 
single-run (no ensemble average over initial conditions) anal- 
ysis for a spatially extended placement of the walker. The pur- 
pose of the latter is to examine the vaUdity of using ensembles 
for single run situations in patch type initial placement. We 
also carry out separately configuration averages of the exact 
and EMT evolutions and find fine agreement. For this latter 
purpose we choose the selfpropagator as the quantity to cal- 
culate. 

We also find that in contrast to the EMT treatment, the exact 
replacement of the disordered system by the ordered system 
with memory, outlined in the introduction and carried out in 
detail in section 5, results in spatially long range memories 
as in earlier analyses of quantum systems Iil2i,il3i1 . This is in 



spite of the fact that the original disordered system has only 
nearest neighbor rates. We display these memories T„ which 
connect a site to another, n sites away. We do this both in 
the Laplace and the time domains, and find that the nearest 
neighbor ones are typically larger by an order of magnitude 
than the others. This explains the success of the EMT even 
though its memories have only nearest neighbor character. 

With the help of our formalism based on Eq. ( fTTI ). we inves- 
tigate finite size effects on effective medium theory and find 
interesting new results: corrections depending on size appear 
in the harmonic mean formula ST4i . A novel result emerges 
involving a bifurcation of the effective long time rate of trans- 
fer Feff as the concentration of broken bonds is varied. The 
bifurcation is transcritical in nature, the vanishing solution for 
Fgff being stable for large concentration of the broken bonds 
and the nonzero solution being stable for small concentrations 
relative to a size-dependent critical value. 

Thus, we have presented a number of extensions of effective 
medium theory in this paper. Not discussed here, but impor- 
tant to point out, are other recent extensions along a line of 
research recently taken by two of the present authors in their 
study of transport on small world networks, particularly of 
the Neumann- Watts kind ll37i [38l [39ll . In those systems stan- 
dard rings (finite chains with periodic boundary conditions) 
with nearest neighbor hopping rates for the random walker 
form the ordered part and additional small world connections 
make up the disordered part. Of particular interest to the de- 
velopments of the present paper is the use of effective medium 
theory to develop memory functions that connect greater than 
nearest-neighbor pairs. Indeed, to correctly describe transport 
on small world networks, as well as on the partially disor- 
dered complex networks of ref. [39], it is generally necessary 
to include memory functions connecting all pairs of sites on 
the network except nearest neighbors, in interesting contrast 
to what we have shown here to be the case for the 1-d dis- 
ordered chain. It is possible that the techniques developed to 
understand complex networks can be applied to understand 
the nature of the spatially long range memory functions for 
disordered systems defined on topologically ordered lattices 
of the sort we have considered in this paper. 

In examining previous work in this field, we find impor- 
tant avenues that were opened by the work of Haus and Kehr 
||23[ I24I [25II . Their approach appears to be similar to ours in 
spirit. Prescriptions exist in their work for going from disor- 
der to explicit forms of the GME or the CTRW. Furthermore, 
their use of projection techniques 1 14] to the disorder problem 
includes important considerations of the initial term L14l [Tsll 
which has been often overlooked ifioll in applications of this 
technique. We have also carefully examined the question of 
the usefulness of the widely quoted analysis of ref. [10]. This 
question is perhaps important given the absence in that anal- 
ysis of a practical prescription for obtaining a usable memory 
function from a quantity describing disorder such as a rate 
distribution function. Our answer is that ref. jlOll helped stop 
the unjustified concerns that some authors [^43] seemed to have 
expressed about the applicability of the GME, equivalently the 
CTRW, to disordered systems. Their message that GME's or 
CTRW's are fully capable of treating disordered systems is 
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correct and valuable. On the other hand, we have explained in 
the Introduction to the present paper, how the correctness and 
applicability of GME/CTRW's can be understood without the 
need for detailed argument. The real need is a practical pre- 
scription for the translation of disorder features into the time 
dependence of memories or pausing time distributions. The 
development and use of such a prescription, already apparent 
in early work 1 16, 24, 25], has been attempted in the present 
paper. 

In concluding, for the use of those who prefer to work 
with CTRW's rather than GME's, we give explicit expressions 
for the CTRW pausing time distribution functions in effective 
medium theory. The situation here (in the EMT, not in the 
exact system as can been seen in Fig. |9]l is separable in time 
and space, and so could be addressed by the formula in ref. 
10]. However we exploit the general relation given first by 
Kenkre and Knox |0,[lil] in Eq. (43) of the first of those refer- 
ences or Eq. (30) of the second. Corresponding to the GIVIE, 
Eq. (|8]l, the CTRW equation which is a sort of non-Markoffian 
Chapman-Kolmogorov equation, becomes iflHl 



with 



P,„(t) = PJO) 



Q„,„(t - s)P„{s)], 
(54) 



Qmni^) — i^m.n+l + ^m,«-l) 

2ne) 



Tie) 



e + ITie) 



(55) 



These formulae can be used after the determination of the 
EMT memory in the various ways we have explained. 
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